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Abstract 

The  electron  stimulated  desorption  ion  angular  distribution  (ESDIAD) 
technique  combined  with  a  digital  subtraction  method  have  been  employed  to 
generate  difference  ESDIAD  patterns  between  measurements  made  at  130K  and 
305K  for  both  adsorbed  Cl  and  OH  on  Si(lOO).  The  difference  ESDIAD  patterns 
reveal,  in  real  space,  images  of  the  thermally  excited  vibrational  amplitudes  of  the 
adsorbates.  This  permits  a  qualitative  visualization  of  the  potential  energy  siuface 
associated  with  an  adsorbate  on  its  adsorption  site.  For  Cl/Si(100),  isotropic 
thermal  broadening  of  the  Cl^-beams  from  130K  to  305K  is  observed.  In  contrast, 
OH/Si(100)  exhibits  am  sotropic  thermal  broadening  of  the  H‘''-beams  on  heating 
from  130K  to  305K.  These  results  are  correlated  with  the  vibrational  frequencies 
of  the  surface  oscillators. 

1.  Introduction 

The  vibrational  dynamics  of  adsorbates  are  an  important  issue  on  surfaces 
since  atomic  motions  are  involved  in  the  elementary  steps  for  surface  adsorption, 
suiface  diffusion  and  surface  reaction  processes.  The  investigation  of  adsorbate 
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vibrational  motions  on  a  well-characterized  single  crystal  surface  can  provide 
detailed  insights  into  the  shape  of  the  adsorbate-substrate  interaction  potential 
well.  The  shape  of  this  potential  well  is  significant  in  governing  adsorption, 
surface  diffusion,  and  surface  reaction.  For  example,  when  a  vibrational  motion 
parallel  to  the  surface  is  thermally  excited  to  a  large  enough  amplitude,  i.  e. 
passing  the  point  of  no-retum  on  the  potential  energy  surface,  the  adsorbate 
vibrational  motion  is  then  co*verted  to  lateral  movement  along  the  reaction 
coordinate  for  either  surface  diffusion  or  surface  reaction. 

ESDIAD  [1,2]  is  a  powerful  technique  for  surface  vibrational  dynamics 
studies,  complementary  to  laser-based  spectroscopy  methods  [3].  The  ESDIAD 
measurement  has  the  advantage  of  directly  imaging  the  thermal  disorder  in  surface 
chemical  bond  directions,  giving  information  about  adsorbate  vibrational  motions 
in  real  space.  In  contrast,  vibrational  spectroscopic  methods  indirectly  reflect  the 
vibrational  movements  in  the  frequency  domain. 

When  an  adsorbate-substrate  system  is  subjected  to  excitation  by  an 
electron  beam,  the  electronically  excited  adsorbate  ejects  molecular  or  atomic 
fi^agments  (ions  or  neutrals)  fi^om  the  surface,  along  a  direction  nearly  parallel  to 
the  direction  of  the  chemical  bond  being  broken.  Since  the  electronic  transition  is 
much  faster  (~10"*^  sec)  than  the  nuclear  motion  (Franck-Condon  transition 
process),  the  spatially  recorded  ESDIAD  pattern  is  a  superposition  of  sub¬ 
picosecond  snapshots  of  the  instantaneous  bond  direction  of  the  adsorbate, 
averaged  over  many  identical  adsorbed  species  undergoing  their  zero-point  and 
thermal  motions.  Therefore,  ESDIAD  provides  a  statistical  picture  of  the 
distribution  of  bond  directions  at  a  particular  surface  temperature. 

Clinton[4]  has  indicated,  from  the  quantum  scattering  analysis  of  the 
Franck-Condon  desorption  cross-section,  that  the  maximum  in  the  ESDIAD 
patterns  corresponds  to  the  average  bond  direction;  the  peak  shape  is  a  reflection 
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of  the  probability  density  function  of  the  components  of  vibrational  motion 
parallel  to  the  surface.  At  a  given  temperature,  T,  the  thermally  averaged 
contributions  of  all  the  nth  vibrational  states,  vj/n(a),  will  give  an  angularly 
dependent  electron  stimulated  desorption  (ESD)  cross  section,  a(a),  at  desorption 
angle  a  as  follows  (with  non-degenerate  vibrational  normal  modes)  [5]; 

a(a)  ~  Z|vi/n(a)|  P(D)  (1) 

a=0 

where  vj/n(a)  is  the  vibrational  wave  function,  and  P(n)  is  the  Boltzmann  statistical 
weighting  factor.  Each  individual  vibrational  state,  contributes  to  the  total 

population  of  the  vibrational  states  at  temperature  T  with  a  statistical  weighting 
factor  of: 


exp  [■(E,-Eo)/lcT] 
2;exp{-(E„-E,)/kT) 


(2) 


Here,  (Eh-Eq)  is  the  energy  difference  between  an  upper  state  and  the  ground  state 
vibrational  levels. 

In  this  report,  we  employ  an  ESDIAD  pattern  analysis  procedure  involving 
the  subtraction  of  a  low  temperature-ESDIAD  pattern  from  a  high  temperature- 
ESDIAD  pattern.  The  difference  pattern  provides  information  about  the  extra 
vibrational  amplitude  excited  at  high  temperature.  A  clear  view  is  presented  in 
real  space  of  the  vibrational  excursions  of  the  adsorbate  in  various  azimuthal 
directions  by  thermal  excitation.  This  information  is  not  available  from 
spectroscopic  methods.  We  show  an  anisotropic  thermal  broadening  of 
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OH/Si(100)  in  contrast  to  an  isotropic  thermal  broadening  of  Cl/Si(100),  consistent 
with  known  structural  and  vibrational  properties  of  the  two  systems. 

The  Si(lOO)  substrate  provides  highly  directional  covalent  bonds  at  its 
surface,  and  is,  therefore,  a  logical  candidate  for  the  study  of  the  isotropy  of 
thermal  broadening  of  surface  bond  directions  by  ESDIAD,  It  is  well  known  that 
the  Si(lOO)  surface  reconstructs  forming  a  (2x1)  structure  from  the  dimerization  of 
the  surface  atoms  [6].  The  reconstructed  surface  consists  of  parallel  rows  of  Si2 
dimers.  Each  surface  Si  atom  has  one  inclined  dangling  bond  pointed  outward 
from  the  surface  which  can  bind  to  the  incoming  adsorbate.  Tlie  sliglit  deviation 
from  the  <100>  direction  after  cutting  a  Si  crystal  results  in  a  real  Si(lOO)  surface 
with  wide  adjacent  terraces  separated  by  monatomic  height  steps,  exhibiting 
orthogonal  Si2  dimer  domains  of  (2x1)  and  (1x2)  structure[7]. 

The  chemisorption  of  CI2  on  Si(lOO)  leads  to  dissociation  [8,9].  At 
saturation,  each  Si2  dimer  is  capped  with  two  Cl  atoms,  one  Cl  on  each  dangling 
bond  inclined  25“  ±  4°  from  the  normal  [8,9].  Previous  studies  have  shown  that 
this  is  the  stable  Cl/Si(100)  structure  produced  by  heating  a  metastable  Cl  layer 
with  a  different  structure  to  673K  [8,9].  The  bonding  structure  and  the  vibrational 
normal  modes  of  Cl/Si(100)  are  shown  in  fig.  I.  Comparing  the  frequencies  of  the 
Si-Cl  in-plane  bending  (180-254  cm’^)  [10,1 1]  and  out-of-plane  bending  (  ~  200 
cm"*)  [11]  modes,  one  expects  a  nearly  isotropic  thermal  broadening  of  the  Cl"^ 
ESDIAD  beams  if  the  bending  vibrations  dominate  the  broadening  process  as 
previously  reported  [12]. 

The  adsoiption  of  H2O  on  Si(lOO)  at  room  temperature  is  also  dissociative 
[13-17].  The  bonding  configuration  for  a  surface  OH  group  is  presented  in  fig.  2, 
together  with  the  approximate  vibrational  normal  mode  frequencies  associated 
with  it  [18-20].  The  large  difference  in  frequency  of  the  SiOH  bending  mode  (~ 
1070  cm'^)  and  the  SiOH  torsional  mode  (~  200  cm'^)  [18-20]  is  consistent  with 
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anisotropic  thermal  broadening  of  the  OH  bond  directions,  since  the  ~  200  cm’* 
mode  will  undergo  extensive  thermal  excitation. 

The  stretching  vibrations  of  Si-Cl  and  0-H  bonds  are  of  minor  importance 
for  the  thermally  induced  broadening  observed  in  ESDIAD  because  of  their 
relatively  high  energy  spacings.  Thermal  excitation  of  these  modes  will  be 
insignificant  at  the  temperatures  employed  here. 


2.  Experimental 

The  ultrahigh  vacuum  (UHV)  system  has  a  base  pressure  of  ~3xl0‘*  *  torr. 
It  is  equipped  with  a  digital  ESDIAD/LEED  apparatus,  a  cylindrical  mirror 
analyzer  (CMA)  Auger  electron  spectrometer  (AES),  an  ion  sputtering  gun,  a 
shielded  quadrupole  mass  spectrometer  (QMS)  for  line-of-sight  temperature 
programmed  desorption  (TPD)  studies,  and  an  additional  QMS  with  an  electron 
gun  and  grid  lenses  for  ion  mass  analysis  during  electron  stimulated  desorption 
(ESD).  A  detailed  description  of  the  UHV  system  and  Si(lOO)  crystal  preparation 
can  be  found  elsewhere  [21,22].  A  multi-capillary-collimated  doser  was  used  for 
CI2  and  H2O  introduction  [23-25]. 

Mass  spectrometric  studies  of  the  ions  produced  during  ESD  indicate  that 
Cl"^  ions  (for  the  Cl/Si(100)  system)  [8,9]  and  H"^  ions  (for  the  OH/Si(100) 
system)  are  the  species  contributing  to  the  observed  ESDIAD  pattern  [26]. 

The  electron  beam  energy  used  for  ESDIAD  studies  is  120  eV.  The  electron 
beam  current,  measured  between  tlie  crystal  and  ground  (uncorrected  for  the 
secondary  electrons),  is  very  small  (~  5  nA),  thus  limiting  beam-induced  surface 
damage  during  data  collection.  Due  to  the  fim'te  detection  solid  angle  of  the 
ESDIAD  apparatus,  a  positive  bias  is  supplied  to  the  crystal  to  compress  the  ion 
trajectories.  This  compression  voltage  is  +10.0  V  for  the  Cl/Si(100)  system  and 
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+100  V  for  the  OH/Si(100)  system.  The  compression  of  ion  trajectories  does  not 
affect  the  qualitative  interpretation  presented  here.  For  all  the  ESDIAD 
measurements,  the  crystal  is  not  heated  by  an  electrical  current,  avoiding  electrical 
effects  on  the  ESDIAD  pattern. 

The  ESDIAD  patterns  at  different  temperatures,  and  after  the  digital 
subtraction,  have  been  azimuthally  smoothed  to  remove  the  slight  inhomogeneities 
of  the  ion  optics  and/or  the  MCP  detector  [27],  The  smoothing  method  is  based  on 
the  symmetry  of  the  observed  ESDIAD  pattern  and  the  assumption  that  ion 
intensities  180®  apart  in  the  azimuthal  angle  should  be  equal  due  to  the  two-fo'd 
symmetry  of  the  Si(100)-(2xl)  and  (1x2)  domains.  To  improve  the  signal-to-noise 
in  the  difference  patterns,  long  data  collection  times  were  employed  involving  the 
collection  of  ~  10^  ions/pattem. 

The  cleanliness  and  structure  of  the  Si  (100)-(2xl)  surface  have  been 
checked  by  AES  and  digital  LEED  for  each  adsorption  cycle.  AES  indicates  a 
clean  surface  with  C(KLL):Si(LMM),  0(KLL);Si(LMM),  and  Cl(LMM);Si(LMM) 
AES  peak-to-peak  intensity  ratios  <lxl0‘3. 

3.  Results 

3.1  Thermal  Broadening  of  ESDIAD  Patterns  For  Cl/Si(100)  and 
OH/Si(100) 

Two  off-normal  Cl"^  ion  beams  are  produced  for  Cl  atoms  bound  to  a 
symmetric  Si2  dimer  site  (fig.  1)  during  ESD  process.  The  probing  of  both  (2x1) 
and  (1x2)  orthogonal  domains  by  the  incident  electron  beam  generates  four  Cl^ 
ion  beams  corresponding  to  four  off-center  peaks  in  the  ESDIAD  pattern.  Figure 
3a  is  a  Cl'^-ESDIAD  pattern  at  I30K.  Figure  3b  shows  an  ESDIAD  pattern 
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obtained  after  warming  up  to  305K  with  all  other  conditions  for  data  collection 
held  constant.  The  broadened  individual  peak  profiles  at  the  elevated  temperatiu'e 
may  be  seen  in  the  perspective  plot  (fig.  3b).  If  the  sample  is  cooled  back  down  to 
130K,  the  low  temperature  ESDIAD  pattern  (fig.  3a)  is  reproduced,  indicating  a 
reversible  thermal  process.  No  chemical  reaction  is  observed  in  this  temperature 
range  [8,9]. 

Reversible  thermal  broadening  of  the  H'*’  ion  beams  is  also  observed  for 
OH/Si(100).  The  OH  group  adsorbed  on  the  dangling  bond  of  a  Si2  dimer  (fig.  2) 
will  produce  an  off-normal  ion  beam  during  the  ESD  process,  and  the  (2x1) 
and  (1x2)  domains  of  Si(lOO)  will  therefore  yield  four  H'*'  ion  beams.  The  130K- 
ESDIAD  pattern  of  ion  beams  from  OH(a)  is  shown  in  fig.  4a  and  the 
broadened  pattern  at  305K  is  shown  in  fig.  4b.  No  chemical  reaction  is 
observed  in  this  temperature  range  [16,17]. 

3.2  The  Difference  ESDIAD  Pattern  For  CI/Si(100) 

While  the  temperature  dependent  ESDIAD  patterns  of  CI/Si(100)  and 
OH/Si(100)  appear  to  be  similar  (fig.  3  and  fig.  4),  tlie  difference  patterns  shown 
in  fig.  5  and  6  are  quite  dissimilar  and  reveal  the  different  dynamical  behavior  of 
Cl/Si(100)  compared  to  OH/Si(100). 

Fig.  5a  results  from  the  subtraction  of  the  nOK-Cl"^  ESDIAD  pattern  (fig. 
3a)  from  the  305K  ESDIAD  pattern  (fig.  3b).  The  total  ion  signal  difference 
pattern  presented  in  fig.  5a  contains  on  the  order  of  10^  €!■*■  counts  for  the  positive 
contribution  matched  by  negative  contribution  of  the  same  magnitude.  The  outside 
boundary  regions  in  the  difference  pattern  are  nearly  zero  in  intensity,  and  provide 
a  secondary  reference  level  for  identifying  the  positive  and  negative  features  of  the 
difference  pattern  (fig.  5a).  The  Cl"''  difTerence  pattern  (fig.  5a)  has  two  features: 
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(1)  The  intensity  gain  at  305K,  corresponding  to  the  positive  intensity  regions  after 
subtraction,  are  found  in  a  strong  normal  peak  at  the  center  of  the  difference 
pattern,  anu  ii  four  smaller  ridges  of  positive  intensity  along  the  diagonal 
directions  (<010>  and  <oTo>  azimuths)  two  of  which  are  centered  on  section  A. 

(2)  The  intensity  losses  at  305K,  corresponding  to  the  negative  intensity  regions  in 
fig.  5a,  are  found  approximately  at  the  Cl^  peak  positions  observed  before 
subtraction  (fig.  3).  Two  of  these  regions  are  centered  on  section  B.  The 
procedures  used  for  generating  the  model  ESDIAD  difference  pattern  shown  in  fig. 
5b  will  be  discussed  later. 


3.3  The  Difference  ESDIAD  Pattern  For  OHySi(lOO) 

In  contrast  to  the  difference  pattern  for  Cl/Si(100)  (fig.  5a),  the  difference 
pattern  for  OH/Si(100)  (fig.  6a)  does  not  show  a  prominent  intensity  gain  in  the 
central  region  at  305K.  Instead,  positive  differences  are  observed  along  ridges  off- 
center  in  diagonal  directions  on  the  contour  plot  (fig.  6a).  Section  A  contains  two 
of  the  four  ridge  regions.  The  ion  signal  presented  in  fig.  6a  has  positive 
contribution  of-  10^  H"*"  ion  counts  and  a  negative  contribution  of  the  same 
magnitude.  The  region  of  negative  contribution  is  located  approximately  at  the  H"*" 
peak  positions  in  the  original  ESDIAD  patterns  (fig.  4).  Two  of  these  negative 
regions  are  located  along  section  B.  The  procedures  used  to  produce  the  model 
ESDIAD  difference  pattern  shown  in  fig.  6b  will  be  described  later. 


10 


4.  Discussion 

4.1  Bond  Vibration  And  Thermal  Broadening 

The  solid  lines  in  fig.  7  give  a  cross  sectional  cut  through  the  difference 
pattern  for  Cl/Si(100)  along  directions  A  and  B  (<010>  and  <01 1>  azimuths)  in 
fig.  5a.  Section  A  in  fig.  7  shows  that  the  central  region  contains  most  of  tlie 
intensity  gain  at  305K.  There  is  also  an  intensity  gain  seen  in  shoulders  on  the 
central  peak.  Since  the  thermally  excited  in-plane  bending  mode  of  the  Si-Cl  bond 
(fig.  lb)  broadens  the  Cl"^  ion  beam,  four  such  broadened  beams  overlap  in  the 
center  region.  In  addition,  the  smaller  intensity  gain  in  the  diagonal  directions 
(<010>  and  <01 0>  azimuths)  of  the  contour  plot  (fig.  5a)  gives  four  off-center 
ridges,  and  two  of  them  are  shown  as  shoulders  in  section  A  (fig.  7).  These  ridges 
are  due  to  the  thennal  excitation  of  the  Si-Cl  out-of-plane  bending  mode  (fig.  Ic), 
which  broadens  the  Cl"^  beam  in  directions  perpendicular  to  the  Cl-Si-Si  plane; 
peak  overlap  fi-om  the  different  Si-Cl  bonds  occurs  here  in  the  diagonal  directions. 

Section  B  displays,  besides  the  intensity  gain  in  the  central  region  of  the 
pattern,  the  intensity  loss  in  the  original  Si-Cl  equilibrium  bond  directions.  The 
thermal  excitation  of  both  the  in-plane  and  out-of-plane  bending  modes  of  the  Si- 
Cl  bond  (fig.  Ib-c)  from  the  ground  state  into  higher  energy  levels  at  305K  leads  to 
the  Si-Cl  bond  having  lower  probability  of  being  found  in  its  equilibrium  position; 
this  yields  a  negative  difference  in  Cl^  intensity  near  this  position. 

The  solid  lines  in  fig.  8  show  the  detailed  cross  section  analysis  of  the 
difference  pattern  for  OH/Si(100)  along  directions  A  and  B  (<0I0>  and  <01 1>)  in 
fig.  6a.  Section  A  shows  that  positive  intensity  is  contributed  primarily  by 
overlapping  features  which  are  located  off  nonnal  along  the  <010>direcfion.  The 
positive  features  in  Section  A  are  due  to  the  thermal  excitation  of  the  Si-OH 


torsional  mode  (fig.  2c)  into  upper  energy  levels.  The  excited  torsional  mode 
broadens  the  H'*’  ion  beam  in  directions  peq)endicular  to  the  H-O-Si  plane.  The 
overlap  of  the  four  broadened  H'^’-beams  fi-om  orthogonal  domains  occurs  in  the 
diagonal  region  producing  ridges  in  the  difference  pattern  (  as  shown  in  the 
perspective  plot  in  fig.6).  Section  B  in  fig.  8  shows  that  the  intensity  loss  at  305K 
is  near  the  original  peak  positions  (fig.  4).  The  same  argument  made  for  the 
intensity  loss  of  Cl/Si(100)  at  305K  applies  for  the  OH/Si(100).  Since  the  thermal 
excitation  of  the  SiOH  torsional  mode  into  upper  levels  depopulates  the  ground 
state,  the  OH  bond  has  a  lower  probability  of  being  found  at  its  equilibrium  bond 
direction,  leading  to  a  reduction  of  H^-beam  intensity  in  this  ducction  at  305K. 
The  origin  of  the  model  curve  will  be  discussed  later. 

4.2  Thermal  Population  of  the  Vibrational  Energy  Levels 

The  qualitative  interpretation  discussed  above  (section  4. 1)  is  supported  by 
the  following  quantitative  analysis.  Using  the  harmonic  oscillator  approximation, 
(En-Eo)  =  nhv,  the  statistical  weigliting  factor  for  the  nth  vibrational  state  in 
equation  (2)  is: 

P(„)=  [exp(-nhv)/kT][l-exp(-hv/kT)]  (3) 

For  the  two  substrate  temperatures  employed  in  the  experiments,  130K  and  305K, 
the  population  for  each  state  at  each  temperature  can  be  calculated  fi^om  equation 
(3)  using  the  estimated  vibrational  frequencies.  Table  1  summarizes  the  population 
of  the  vibrational  state  for  the  Si-Cl  and  Si-OH  normal  modes.  The  stretching 
vibrations  for  the  Si-Cl  and  the  0-H  bonds  are  neglected  due  to  their  minor 
relevance  to  the  thennal  broadening. 
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For  the  Cl/Si(100)  system,  the  following  conclusions  are  evident  from 
analyzing  the  population  of  each  vibrational  state  as  a  function  of  the  substrate 
temperature;  (1)  There  is  a  substantial  depopulation  (AP(n=0)  ~  28%)  of  the 
ground  state  for  both  in-plane  and  out-of-plane  bending  modes  of  the  Si-Cl  bond  at 
305K  compared  to  130K.  The  depopulation  will  lead  to  an  intensity  loss  upon 
heating  in  the  equilibrium  Si-Cl  bond  direction  as  observed  in  the  difference 
pattern  (fig.  5a);  (2)  At  305K,  the  n  ^  1  vibrational  states  for  both  the  in-plane  and 
out-of-plane  Si-Cl  bending  modes  are  considerably  populated  as  shown  in  Table  1. 
The  population  of  upper  states  will  result  in  a  larger  vibrational  amplitude  in  the 
directions  associated  with  these  modes,  as  discussed  later  (section  4.3);  (3) 
Isotropic  broading  of  the  Si-Cl  pattern  is  expected  based  on  the  presence  of  in¬ 
plane  and  out-of-plane  vibrational  modes  of  similar  frequency  (~  20^  cm"^). 

Different  vibrational  dynamics  are  expected  for  the  SiOH  bond.  (1)  The 
large  energy  for  excitation  of  the  SiOH  bending  mode,  ~  1070  cm‘^  [18-20],  leads 
to  little  excitation  of  the  SiOH  bending  mode  over  the  130K  to  305K  temperature 
range.  (Table  1).  Therefore,  the  broadening  in  the  azimuthal  direction 
corresponding  to  the  H-O-Si  plane  (fig.2)  is  expected  to  be  small  for  the 
temperature  range  studied  here  (130K  -  305K);  (2)  The  SiOH  torsional  mode  will 
be  the  major  contributor  to  shape  changes  in  the  H'^  ESDIAD  pattern  on  raising 
the  temperature  because  of  its  low  frequency;  (3)  The  factors  mentioned  in  (1)  and 
(2)  will  result  in  an  anisotropic  thermal  broading  effect. 

4.3  Mean  Square  Amplitude  of  Normal  Modes 

The  mean  square  amplitude  of  an  oscillator  is  directly  related  to  its  quantum 
state  specified  by  the  quantum  number,  n,  within  the  hannonic  oscillator 
approximation.  Using  the  normal  coordinate,  Q  =  m’^^  x,  as  discussed  by  Wilson, 
Decius  and  Cross  [28],  with  x  as  the  displacement  and  m  as  the  mass,  the  mean 
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value  of  an  observable  (the  square  amplitude,  <Q^>,  in  this  case)  in  its  quantum 
state,  can  be  evaluated  as  129,30]; 

<QS  >  =  (H^n(Q)|Q^|v«(Q))=  I  ^ 

in  which  u  is  the  wavenumber  of  the  normal  mode  and  c  is  the  speed  of  light. 
From  eqn.  (4),  the  mean  square  amplitude  for  a  given  normal  mode  is  expected  to 
increase  linearly  with  the  quantum  number,  n.  Furthermore,  it  is  evident,  from 
eqn.  (4),  that  a  low  frequency  oscillator  will  have  a  larger  mean  square  amplitude 
than  a  high  frequency  one. 

4.4  Modeling  of  The  ESDIAD  Patterns  For  Cl/Si(100)  and  OH/Si(100) 

Both  Cl/Si(100)  and  OH/Si(100)  systems  have  been  modeled  for  the 
vibrational  behavior  of  the  Si-Cl  and  SiOH  bonds.  In  the  modeling  of  Cl/Si(100), 
the  following  optimized  parameters  have  been  employed  for  the  305K  ESDIAD 
pattern  compared  to  the  130K  ESDIAD  pattern  using  a  Gaussian  distribution  for 
the  individual  ion  beams  profiles:  (1)  an  isotropic  20%  increase  in  full  width  at 
half  maximum  (FWHM)  of  the  Gaussian  distribution;  (2)  a  10%  increase  in  the 
integrated  ion  yield  corresponding  to  that  measured  experimentally;  (3)  a  5% 
decrease  in  the  apparent  Cl"^  ion  polar  angle  detected,  which  is  not  explored 
further  here.  The  modeled  difference  pattern  (fig.  5b;  fig.  7  with  dashed  lines) 
agrees  quite  well  with  the  experimental  results  (fig.  5a;  fig.  7  with  solid  lines). 

For  the  OH/Si(100)  system,  the  following  parameters  are  used  for  modeling 
the  ESDIAD  pattern  at  305K  compared  to  130K:  (1)  a  2%  increase  of  the  FWHM 
of  the  beam  Gaussian  distribution  in  the  Si2  dimer  bond  azimuthal  direction 
component  (in  0-Si-Si  plane).  In  contrast,  an  8%  increase  of  the  FWHM  is 
employed  for  the  orthogonal  azimuthal  direction  component;  (2)  an  increase  in  the 
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integrated  volume  of  the  Gaussian  peak  profiles  of  -3%,  in  accordance  with  the 
total  H"*"  ion  yield  measured  experimentally;  and  (3)  a  small  decrease  (~6%)  in 
polar  angle  for  the  beams,  for  which  the  origin  is  not  explored  here.  The 
modeled  difference  pattern  (fig.  6b;  fig.  8  with  dashed  lines)  also  agrees  well  with 
the  experimental  results  (fig.  6a;  fig.  8  with  solid  lines).  Thus  in  this  case  the 
anisotropy  for  the  two  -OH  modes  is  approximately  a  factor  of  four. 

It  is  not  appropriate  in  the  analysis  to  quantitatively  estimate  (Qu)  fron* 
ESDIAD  data  presented  here,  since  the  laboratory  angles  measured  depend  upon 
crystal  bias  and  ion  energy  and  therefore  possess  a  complex  relationship  to  (On)- 

The  above  analysis  of  the  vibrational  dynamics  of  the  adsorbates  has 
qualitatively  explained  the  observed  isotropic  and  anisotropic  thermal  broadening 
processes.  However,  it  should  be  kept  in  mind  that  there  are  additional  factors  to 
be  evaluated  for  a  full  interpretation  of  the  results,  such  as  the  beam  broadem’ng 
caused  by  the  retention  of  the  momentum  of  the  initial  state  by  the  departing  ion  of 
the  adsorbates  [5,  3 1],  and  the  perturbations  of  the  departing  ion  motions  (final 
state  effect).  Two  prominent  final  state  effects  are:  (1)  Ion  trajectory  distortion 
toward  the  surface  by  an  electrostatic  image  potential  (a  long  range  effect);  and  (2) 
Ion  trajectory  distortion  away  from  the  surface  due  to  preferential  neutralization  of 
the  desorbing  ions  in  the  near  surface  region  (a  short  range  effect).  These  two 
effects  have  been  studied  in  detail  elsewhere  [8,  32-33]. 

5.  Conclusion 

The  ESDIAD  method  is  shown  to  be  a  liigh  resolution  technique  for  study 
of  adsorbate  vibrational  dynamics.  By  employing  ESDIAD  pattern  subtraction 
methods,  the  tliermal  disorder  in  chemical  bond  directions  on  the  surface  may  be 
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mapped,  giving  a  qualitative  insist  into  of  the  potential  energy  surface  for  an 
adsorbate  near  its  equilibrium  bonding  position.  Examples  of  isotropic  vibration 
(Cl/Si(100))  and  anisotropic  vibration  (OH/Si(100))  are  presented  and  the  relation 
to  the  estimated  normal  mode  frequencies  has  been  established. 
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Figure  Captions 


Fig.  1.  Configuration  and  normal  modes  of  terminally  bonded  Cl  atoms  adsorbed 
on  a  symmetric  Si2  dimer  site  on  Si(lOO). 

Fig.  2.  Surface  species  resulting  from  H2O  dissociative  adsorption  on 

Si(lOO)  with  the  OH  on  one  dangling  bond  of  the  symmetric  Si2  dimer 
and  the  H  atom  on  the  other. 

Fig.  3.  Reversible  temperature  dependent  ESDIAD  patterns  of  Cl/Si(100).  (a) 
130K  Si  substrate  temperature;  (b)  305  K  Si  substrate  temperature.  The 
figiues  on  the  left  are  perspective  plots;  they  are  generated  by  mapping  the 
counts  of  Cl"^  ions  in  the  z  direction  as  a  function  of  the  (x,y)  coordinates 
on  the  detector.  The  figures  on  the  right  are  contour  plots.  Each  contour 
line  in  each  individual  plot  represents  an  increment  of  1/7  of  the  maximum 
count  rate  at  a  peak,  with  a  width  of  ±5%  of  the  median  value  of  that 
contour. 

Fig.  4.  Reversible  temperature  dependent  ESDIAD  patterns  of  OH/Si(100).  (a) 
130K  Si  substrate  temperature;  (b)  305  K  Si  substrate  temperature.  H2O 
was  adsorbed  at  305K  to  saturation.  The  figures  on  the  left  are  perspective 
plots;  they  are  generated  by  mapping  the  counts  of  H"*"  ion  of  the  adsorbed 
OH  group  in  the  z  direction  as  a  function  of  the  (x,y)  coordinates  on  the 
detector.  The  figures  on  the  right  are  contour  plots.  Each  contour  line  in 
each  individual  plot  represents  an  increment  of  1/7  of  the  maximum  count 
rate  at  a  peak,  with  a  width  of  ±5%  of  the  median  value  of  that  contour. 
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Fig.  5.  The  difference  ESDIAD  patterns  for  Cl/Si(100).  (a)  Experimental 
difference  pattern  obtained  by  digital  subtraction  of  130K  ESDIAD 
pattern  (fig.  3a)  from  305  K  ESDIAD  pattern  (fig.  3b).  (b)  Modeled 
difference  ESDIAD  pattern  with  Gaussian  distribution  peak  profiles. 

Fig.  6.  The  difference  ESDIAD  patterns  for  OH/Si(100).  (a)  Experimental 
difference  pattern  obtained  by  digital  subtraction  of  130K  ESDIAD 
pattern  (fig.  4a)  from  305  K  ESDIAD  pattern  (fig.  4b).  (b)  Modeled 
difference  ESDIAD  pattern  with  Gaussian  distribution  peak  profiles. 

Fig.  7.  Cross  section  views  of  the  difference  ESDIAD  pattern  of  Cl/Si(100) 
shown  in  fig.  5  cut  in  directions  A  and  B.  Solid  lines  represent  the 
experimental  data  while  the  dashed  lines  show  model  calculations 
involving  isotropic  broadening  of  the  four  Cl"''  ESDIAD  beams.  The 
angles  presented  in  the  figure  are  the  measured  angles  without  correction 
for  ion  trajectoiy  and  final  state  effects. 

Fig.  8.  Cross  section  views  of  the  difference  ESDIAD  pattern  of  OH/Si(100) 
shown  in  fig.  6  cut  in  directions  A  and  B.  Solid  lines  represent  the 
experimental  data  while  the  dashed  lines  show  model  calculations 
involving  anisotropic  broadening  of  the  four  H"''  ESDIAD  beams.  The 
angles  presented  in  the  figure  are  the  apparent  angles  without  correction 
for  ion  trajectory  and  final  state  effects. 


Schematic  of  The  Normal  Modes 
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Schematic  of  The  Normal  Modes 
of  OH  on  Si{ioo)-(2xi) 
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ESDIAD  Patterns  of  Cl  on  Si(lOO) 
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ESDIAD  Patterns  of  OH  on  Si(lOO) 
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Difference  ESDIAD  Patterns  of  CySi(lOO) 

(a)  Experimental 


Difference  ESDIAD  Patterns  of  OH/Si(100) 
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(b)  Model 
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Cross  Section  Views  of  The  Difference 
ESDIAD  Patterns  For  Cl  on  Si(lOO) 
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